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presented also for this theory. It seems to support, as well, physical non-equivalence of the 
two frames. 
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I. INTRODUCTION 



One of the most important problems of modern cosmology (and theoretical physics too) is the 
explanation of the current universe speed-up, first discovered in A convenient way to express 

this situation is to introduce a new form of energy, called Dark Energy (DE). Indeed, it is not 
easy to generate the necessary amount of repulsive force (not less the very fact that it is repulsive) 
through quantum vacuum fluctuation contributions of ordinary fields, at cosmological scale f3]. To 
start with, one needs to find an acceptable solution to the cosmological constant problem. During 
the last few years several theories, which are extensions of or alternatives to Einstenian Gravity, 
have been developed in order to formulate and try to explain the dark energy universe. The most 
accurate observational data we now have indicate that the equation of state (EoS) parameter, u, 
for DE is very close to —1 (for a review of observational data from the theoretical point of view, 

n n 

see e.g. [4J, and for a description of the observable cosmological parameters, see e.g. [5(]). 

As gravitational alternatives for DE, modified gravity theories have been formulated, calling 
for plausible late-time modification of General Relativity (GR). Many modified gravity models 
have been proposed in the literature (for a review, see [g]), starting from the very simple 1/R 



theory [7 



(that soon was declared as problematic) to more elaborated ones which, being still 
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not fundamental, are already quite often inspired by string and M-theory considerations 

There are some modified gravity theories (see e.g. Q]) which have the capacity to achieve 
a reasonably natural unification of the early-time inflation epoch with the late-time acceleration 
period, in which we are now, by taking advantage of the different role played by the gravitational 
terms relevant at small and at large curvatures. Moreover, if one assumes that the universe may 
enter in a phantom phase (lo < —1) [12II31I]. modified gravity may naturally describe the transition 
from the ordinary, non-phantom phase towards the phantom one without the necessity to introduce 
exotic matter. Also, no Big Rip is usually expected there, because in modified gravity the phantom 
phase is often transient. These theories may also be able to encompass, in a truly unified way, DE 
and dark matter. 

For any such theory to be valid it is always strictly required that it accurately describes the 
known sequence of cosmological epochs, specifically it must fit very well an increasing number of 
more and more precise observational data fl^l fis|. In particular, it must pass the Solar System 
tests (for f(R) models, this is discussed e.g. in Refs. [21, 



while the case of SdS metrics is 



studied in [19|). After all this procedure, it turns out, in the end, that the possibility exists to pass 
all known local tests and to construct a realistic f(R) modified gravity theory which unifies 
early-time inflation with late-time acceleration [13I ]. As a consequence, a natural resolution of the 
problem at hand seems to be on the right track. 

Another important issue to be taken into account concerns the appearance of different types of 
finite-time, future cosmological singularities (see [21]), what is a typical property of the number of 
modified gravities with the effective quintessence/phantom behavior, in the same way as it had been 
found for other, more simple scalar/fluid dark energy models. It has been proven (see jzl [23]) 
that in modified gravities all four known types of finite-time future singularities may possibly 



appear. However, modified gravity has the chance to cure all of these future singularities b y th e 
addition of an R 2 -term. which is relevant at the very early stages of the universe only (see [241]). 
The combination of non-singular realistic modified gravity with singular realistic models unifying 
the early-time inflation epoch with late-time acceleration finally results in a realistic, non-singular 



theory (for more details see 25(|). For all these reasons, modified gravity theories constitute a very 
attractive alternative to dark energy. 

In this paper, by making one further step in the direction of trying to build a truly realistic 
theory, we develop a non-trivial variant of the accelerating cosmology reconstruction program both 
for f(R) gravity and for a modified Yang-Mills/Maxwell theory (see [27]] for related work). In 
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Sect. II, two seemingly different reconstruction schemes, (i) in terms of e-foldings (for a general 



review, see 



2a]), and (ii) by using an auxiliary scalar field (see, e.g. 28B31I]) are reviewed and then 
explicitly compared, what is done here for the first time. To illustrate the results, the example 
of a model with a transient phantom behavior without real matter is discussed in both schemes. 
In Sect. Ill the corresponding reconstruction schemes for the case of a Yang-Mills theory (in a 
extended, non-minimal version) are developed, one of them using an auxiliary scalar field and the 
other one without it. The specific example of power-law expansion is carefully considered in both 
schemes. Finally, in the last section, a summary of the results obtained and, in particular, of the 
detailed comparison with pros and cons of the two schemes of reconstruction is presented. 

II. COSMOLOGICAL RECONSTRUCTION OF REALISTIC MODIFIED f(R) 

GRAVITIES 

A. GENERAL FORMULATION IN TERMS OF E-FOLDINGS 

We will shaw in this section how one can construct an f(R) model realizing a given cosmology, 



by using the techniques of [32|. For the benefit of the reader (and self-consistency of the paper), 
in a simpler situation, and with the help of an explicit example, we shall review the method that 
will be later applied to the Yang- Mills case. The starting action for f(R) gravity (see e.g. jfj], for 
a general review) is 



S 



/ d X^J-g I + ^matter I • (1) 



The first FRW equation turns into the following field equation 

= - + 3 (H 2 + H) f'{R) - 18 Uh 2 H + HH) f"{R) + K 2 p, (2) 



2 

with R = 6H + 12H 2 . Using the e- folding variable, N = In instead of the cosmological time t, 
one gets 

= -1^1 +3(H 2 + HH') f'(R) - 18 (4H 3 H' + H 2 {H') 2 + H 3 H") f"{R) + n 2 p, (3) 

where H' = Assuming the matter density p is given in terms of a sum of fluid densities with 
constant EoS parameters, w,, we have 

P = £ P m^) = ^ p l0 ao 3(1+ ^e- 3{1+ ^ N . (4) 
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Using the Hubble rate H = g(N) = g{— ln(l + z)), with z = e — 1 the redshift, the scalar 
curvature takes the form: i? = Qg' (N)g(N) + 12<7(A^) 2 , which can be solved with respect to N as 
N = N(R). Defining G(N) = g(N) 2 = H 2 and using g]), Eq. © yields 



= -9G(N(R)) [4G'(N(R)) + G"(AT(i2))] + 

+ ^ /Oj0 G - 3(1+ ^ ) e- 3 ( 1+ ^) 7V ^. (5) 

i 

This is a differential equation for f(R), where the scalar curvature is here R = 3G'(N) + 12G(N). 

1. EXAMPLE: ASYMPTOTICALLY TRANSIENT PHANTOM BEHAVIOR 
Let us consider an evolution given by the following Hubble parameter: 



H\N) = H In 1-1 +H t = H N + H X = G(N), (6) 

where Hq and H% are positive constants. We can, in this case, achieve a phantom behavior with 
the possibility to be asymptotically transient, without the presence of real matter. The present 
is actually a simplified example, but can be of use as a component part of a more elaborated 
model where, with a modified functionality, the transition can be reached at finite time (this is 
work in progress, the results of which will be presented in a future publication.) Indeed, from 
R = 3G'(N) + 12G(iV), we find 

12Hq Hq 

Eq. (jll Ap takes the form (with r being the curvature measured in terms of Hq, r = R/Hq) 

and changing now the variable from r to x, as x = ^jsp, Eq. ([8]) reduces to 

d 2 F(x) ( 1\ dFix) 7— i / \ 

which is easily recognized as a degenerate hyper geometric equation, whose solutions are given by 
the Kummer's series 3>(o, b;z), the simplest one being 

1 r-3\ „ / 1 r 2 



f{r) = C*{-2,--;— (10) 

where C is a constant. As a consequence, with this f(R) theory, the solution given by Eq. (I10p . 
we can reproduce the phantom behavior without real matter given by Eq. (|6|). 
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Taking this into account, for (|6|), we have 

H(t) = ^(t-t ), (11) 

and it turns out that, with this model, we have a contribution of an effective cosmological constant 
and another term which will produce an accelerating phase but, remarkably, without developing a 
future singularity, in spite of its phantom nature. Hence, the f(R) gravity given by Eq. fjlOj) gives 
rise to a cosmological solution, with an asymptotically transient phantom behavior, which does not 
evolve into a future singularity. This property relies on the fact that the phantom behavior gets 
more and more mild with time (asymptotically disappears), at a rate that overcomes the one for 
the formation of the singularity. 

Actually, there is another independent solution of Kummer's equation ([9]), the complete solution 
being: 



where the second basic solution is a Laguerre L function, which is well behaved but cannot be 
represented as a simple rational one. It is interesting to note that this second function asymptoti- 
cally behaves exactly in the same way as the first, for large negative curvature (e.g., as R 2 , when 
R —7- — oo). For large positive one it explodes exponentially, as R~^l 2 e R / 12 (again, R in units of 
Ho). 

B. GENERAL FORMULATION USING A SCALAR FIELD 

In this section it will be shown how to construct an f{R) model realizing a given cosmology, 
but using this time a different technique, which involves a scalar field [331 ] . The final aim will be to 
apply this procedure to the novel case with a Yang-Mills term, what will be performed in Sect. III. 
Here we summarize the basic tools necessary in order to understand the procedure and to make 
the present paper self-contained. We start from the action for modified gravity 



S = J d 4 x^(f(R) + C matter ), (13) 

which is equivalent to 

S = J d A x^{P(ct>)R + Q(<f>) + C matter ). (14) 

Here, C ma tter is the matter Lagrangian density and P and Q are proper functions of the scalar 
field, (j), which can be regarded as an auxiliary field, because there is no kinetic term depending on 
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(j) in the Lagrangian. By varying the action with respect to </>, = P'{(p)R + Q'(cp), which can be 
solved in terms of (f>, as <j> = 4>(R)- Substituting it into (I14p and comparing with (|13h . one obtains 

S = j ' d?Xy/=jj(f(R) + Crnatter), f(R) = P((p(R))R + Q((p{R)), (15) 

and varying the action with respect to the metric g^, 

= ~g^{P{<f>)R + Q{4>)) + R^P(4>) + g^V 2 P(^) - V M V„P(0) - \t^. (16) 
The equations corresponding to the standard, spatially-flat FRW universe are 

= -Q(4>) - 6P 2 P(0) - 6P^^ + p, (17) 

= Q(J>) + (AH + 6P 2 ) P(0) + 4F^M + + p, (18) 

and, by combining them, we find 

dt A at 

As we are allowed to redefine the scalar field <f> properly, we choose the most simple, non-constant, 
smooth possibility (what is commonly done in this kind of approaches), namely <fi = t. 

Now, given a cosmology, specified through the scale factor a, given by a proper function g(t) as 

a = a e 9(t \ (20) 

with a constant do, and if it is assumed that p and p consist of the sum of different matter 
contributions, each one with constant EoS parameter, u)i, then Eq. f)19|) reduces to the following 
second order differential equation 

= 2 ^ " 2 ^^T + + EC 1 + ^)PiO% 3{1+U}i) e-^ + ^\ (21) 



from where one can obtain P(<j>) and, using Eq. (fT7|) . 

Q(0 = -6( 5 '(^)) 2 P(^) - 6<M)^ +J2pi0% 3il+ ^e-^+^. (22) 

i 

As a result, and as anticipated, any given cosmology (|20p can indeed be realized through some 
corresponding /(i?) -gravity. Before applying this procedure to the Yang-Mills case, let us make 
things even more clear by means of the example considered before. 
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1. EXAMPLE: ASYMPTOTICALLY TRANSIENT PHANTOM BEHAVIOR 



In order to compare the two different methods developed for the reconstruction of f(R) 
gravities — to reproduce any given cosmology — we consider again the asymptotically transient phan- 
tom behavior, without real matter, given by ([6]). The Hubble parameter can be written as 



H = yjH g{t) + Hy 



dg(t) 
dt 



and thus 



H 



Hi 



with c an integration constant. Introducing (|24p into (|2ip . 

d 2 P(cft) H ^P(<f>) + rr pU , 

and using a new variable, x = (ft — c, we get 

d 2 P(x) H Q dP(x) 







whose solution is 



dx 2 2 dx 
( 



+ H P(x), 



(23) 



(24) 



(25) 



(26) 



P(x) = \ (2 - H x 2 ) d + \ (2 - Hox 2 ) C 2 



\ 



e 4 x x 



4(2 - # x 2 ) 4^ 



e y dy 



\ 



(27) 



/ 



Now, using ([IT]), we obtain 



Q(x) = ^ H o x 



8H x (2 + Hox 2 ) Ci 



2 /* 

(8 + 2^ ^ 2 ) + i 2y / #ox (2 + i^ 2 ) / e~ y2 dy 

o 

(28)' 



Co 



Taking into account that R = 6H + 12# 2 = 65" (x) + 12 (g'(x)) , it follows that 



'R-3H 
3H 2 ' 



(29) 



Introducing at this step (|27p and (|28|) into (|15p . and considering (|29p . one finally gets the explicit 
expression 



f(R) 



R 2 - 18H R + 9# 2 
V2Ho 



12H 



16 24H 







0; / e y dy 





(30) 



C 2 



We thus have proven that, within this scheme, we are able to obtain the f(R) model (|30p which 
reproduces the desired transient phantom behavior without real matter, as given by ([6]). 
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III. RECONSTRUCTION OF THE YANG-MILLS THEORY 



GENERAL FORMULATION 



In this section we develop a reconstruction scheme of the YM theory without any auxiliary 
scalar field. Consider the following action 

R 



S 



dx q 



-9 



2 K 2 ~ \ V ) 



(31) 



where F® = d^Af, — d u A°^ + f abc A b ^A1, and T may be assumed to be a smooth function (however, 
strictly speaking it may suffice if it is continuously differentiable) . The presence of this function 
is necessary in order to allow for more freedom in the choice of the theory 31[ (since difficulties 



inherent to the problem may prevent obtaining the standard Yang-Mills case). Anyway, this 
function will not constitute a problem for the development of our methods, which are thus proven 
to be even more powerful. For simplicity of the derivation (and in order not to break the line of 



argument of this paper) we concentrate here on the SU(2) case where f abc 



r abc 



but, with some 



more effort, exactly the same procedure can be extended to other gauge groups (more general cases 
will be treated in a subsequent publication). Taking into account that 



S{F^F a ^) 



S ( d a A\ 



ha/3 



4F 



the equation of motion for the field potential A a turns into 



(32) 
(33) 



5S 



5 {d v A°) 



and, from here, 



-gT'iF^F^ 



pa v^i 



+ 



ss 



-gF'{F^F a ^ 



^abc j^b pc ufj, 



0. 



(34) 



(35) 



Variation of (|3ip with respect to g^ v yields the following equation of motion 

± (iV - \g^R} - {ifyF"*) + 2F' {tyF"*) F^Fy = o, (36) 

S( F a F a P a \ / \ 

where we have used v ; = 2F' [F^F aa ^\ F^F? 1 . Considering now a FRW universe, and 



the following Ansatz for the gauge field 



3J 



A a 



H = 0, 



(37) 
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the /i = component of f|34[) becomes an identity, the ji = i component yields 



d t \a(t)F' (F^F aa f ) ) A(i)e A W 
while the (t, t) component of ([36]) is 

3ff2 ( t ) , I jr ( 77a P ac>/3 
2k 2 ^ 2 V 

and the (i, i) component of (f36|) reduces to 



+ jr' ( pa F aaf}\ 3A(t) = q 

a(i) v ^ ; 



/ \ X 2 OV 2A (*) 

K? FaaP ) + 2 <* 2 F' [F^F aa ^ ) \^ = 0, 



1 

2^ 



2ij(t) + 3# 2 (f) 
Adding ([39]) to (|30"]l . one arrives at 



2A(t) 

2 ' ' T'(F^ R F aat 



a 2 (t) 



X 2 (t) - 2a 2 



(38) 



(39) 



e 2A(t) 



0. (40) 



(t)H(t) - 4 K 2 a 4 F' (F^F aa ?) e 4A W = 0, 



and then 



F' {F^F 11 '' : ) - ° ( f ) c -4Affl 



4/€ 2 a 4 



(41) 



(42) 



Using flU, Eq. fl38|) reduces to: 



2a 2 H(t)e 2X ^ 



5a(t)a(t)H(t) + a 2 (t)H(t) X(t) - 3a 2 (t)H(t)X 2 (t) + a 2 (t)H(t)X(t) = 0, (43) 



which constitutes a differential equation for X(t). Hence, by using Eq. (|37p . once we have the 
function X(t), given by (|43p . we can obtain the corresponding YM theory which reproduces the 
selected cosmology. The Ansatz considered above actually leads to a mathematical solution of the 
problem. 



1. EXAMPLE: POWER LAW EXPANSION 



t ^ 



Considering the case of power law expansion: a(t) = y^j > where t\ and h\ are constant, 
and assuming X(t) = (hi - I) In (±) + Ai, where Ai is again a constant, Eq. (|43p reduces to the 



following algebraic equation 

hi(hi - 1) + a 2 t 2 e 2Xl = 0, (44) 

hence 
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and 



A(*) = (*i-Dm(l) + iiB^-^ 



aH 2 



(46) 



With the help of this reconstruction scheme, we have obtained the function X(t) given by ([46 p . 
Then, using (|37p . we are able to reproduce the cosmology given by the power law expansion: 
a(t) = 



B. GENERAL FORMULATION USING A SCALAR FIELD 

In this section, the reconstruction of a YM theory is performed using the technique of [3 1 . By 
way of introducing an auxiliary scalar field, eft, we can rewrite action (|31|) as 



(47) 



Variation of (j47]) with respect to <j> yields the corresponding equation of motion 

n pa pa^dPtt) dQ{4>) 

which can be solved with respect to <f> as <f> = 4>(F"' lJ F a ^' / ). Then, action (|47p is rewritten as 



(48) 



S 



dx \J—g 



R 



Ik 2 ^ ^ v > 



(49) 



where 



jr^pa^ = l P ^ F a vF a n)F a vF a^ + In^pa^pa^ 



(50) 



^ \T V jW II ilf ' £ 

Taking the variations of this action (I47j) with respect to g pu , we obtain the Einstein equation 

^2 ( R v» ~ l^vRj = -\p{^)F a , p Fr + l 9lw (pWlfyF"* + Q{4>)) . (51) 

Finally, taking the variations of ([37]) with respect to A^, it follows that 

= d u (V=gP(^)F au ») + ^P{4>)f hc A h v F CVi " . (52) 

We restrict our analysis to the case where the gauge algebra is SU (2) and the gauge fields are given 
by 



A a 



0, 



n = o. 



(53) 
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With these assumptions, Eq. (jiS|) reduces to 

= 6 (-a 2 X(t) 2 e 2X ^a(t)- 2 + a 4 e 4X ^a(t)-^ dP(0) 



+ 



The (t, t) component of (f51~|) is 

= + I (a 2 \{tfe 2 ^a(t)- 2 + a 4 e 4A «a(r 4 ) Pfo) + ±Q(0), 

k z 2 V / 4 

the (i, i) component of (|5ip becomes an identity, while the component is 



(54) 



(55) 



2k 2 



2H(t) + 3if(t) : 



P{<t>) fa2^ (t) 2 e 2A(t)«(t)-»+«^Wa(t)- 



+ 



<5«. (56) 



The /U = component of (|52p becomes an identity, and the \i = i component yields 

= d t (a(i)P(</>)A(t)e A W) + 2a 2 a(t)- 1 P(cf ) )e 3X ^ . (57) 

Here, we can identify <f> = t, because we are always allowed to take the scalar field <f> properly in 
order to satisfy this. By differentiating (|55p with respect to t and eliminating Q = , it follows 
that 

= 2-H(t)H(t) + a 2 X(t) 2 e 2X(t) a{t)- 2 P(t)+ 



+ 



a 2 (A(i)A(i) + X(t) 3 - X(t) 2 H(t)^j e 2X ^a{t)~ 2 + 2a 4 (a(*) - H(tj\ e 4X ^a(t)- 4 ^ P(t). (58) 



Using ([57]) . we can solve for P in (jUIBp . and obtain 

a(t) 2 H(t) 



K 2 a 2g2A(t) 



A(i) 2 + a 2 e 2A Wa(t)- 2 



(59) 



Taking into account (|59j) . Eq. (|57l) reduces to 



= 2H(t) (a 2 e 2X ^a(t)- 2 ) + a 2 e 2X ^ a(t)~ 2 \\(t) (&H(t)H(t) + H(tj) - X(t) 2 H(t) + X(t)H(t) 



+X{t) 2 X(t) 3H(t)H(t) + H(t) - X(t) 2 H(t) - X(t)H(t) 



+ 



(60) 

which constitutes a differential equation for X(t). As was the case for the other reconstruction 
scheme, developed in section III. A, once we have the function X(t) — given here by QIIIBp — we can 
readily obtain the modified YM theory that reproduces the desired cosmology, through the use of 
Eq. (j53p . which was our starting Ansatz. Note that in Oil Bp we have positively corrected some 
missprints of a previous calculation (recognized by the authors). Using Eq. (|55j) and (f59j) . it is 
easy to check that Eq. ([56]) becomes an identity, thus being always fulfilled. 
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1. EXAMPLE: POWER LAW EXPANSION 



Consider now the case of power law expansion: a(t) = ijr) > where t\ and h\ are constants 
and assume X(t) = (hi — 1) In ( j- J + Ai, where Ai is a constant. Eq. Oil Bp turns into the algebraic 



one 



hence 



and then 



= aHfe 4Xl + a 2 t?e 2Al (/ii - l)(2/ii - 1) + %(^i - l) 3 , 



(61) 



Ai = < 



I In f Mi^M 



L 2 



272— , 



(62) 



A(i) = [hi - 1) In ( - ] + < 



i ln ( fei(i-fei) 



k 2 



(63) 



With the function X(t) given by (|6^|) and using ([53^1 . we can now reproduce the cosmology given 
by the power law expansion: a(t) = 



IV. SUMMARY AND DISCUSSION 



In this paper, after carefully reviewing and comparing (for the first time), and also with the 
help of corresponding examples, two different schemes of reconstructing cosmologies for modified 
gravity, we have successfully extended them to the case of YM theories. The first scheme does 
not need an auxiliary scalar field, while the second one is thoroughly based on its use. With 
these reconstruction methods, any explicitly given cosmology can be realized as a corresponding 
modified gravity or YM theory, respectively. Although this fact had been already anticipated in 
the specialized literature — for the first of the two situations and concerning some basic models — it 
is comforting to see here how it can be also explicitly extended to more realistic physical theories, 
as the modified Yang-Mills one, with reasonable effort. As we have indicated in the paper, things 
are far from straightforward, and a very careful analysis of the solutions here obtained (also in 
relation with the comparison of the two different methods) , and of other additional solutions of the 
differential equations, with potential physical interest, is still necessary. First results indicate that 
the solutions obtained do pass the solar system tests and the other known physical constraints. 
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In order to compare both schemes of reconstruction, an example has been explicitly worked out 
in the two cases. The result obtained in the first scheme (in terms of e-foldings) is 



f(R) = C$ 



1 R-3H 
' ' 2' 12H 



C 



1 1 

~4 + 2#o 



R 



1 



36# 2 



R* 



(64) 



while for the second scheme (using an auxiliary scalar field), the f{R) obtained has the form 

(R 2 - 18H R + 9H$) 



f(R) = - R2 - l8 "f + 9H °c 1 



l2H n 



iL_Ma R - 9H 



R-3H 
12H 



16 



24# n 



-y 



dy 



(65) 

As one can easily see, the results obtained for both methods are in fact different. The reason behind 
this is the fact that action (|14p corresponds to a wider class of theories than action (jl3j) (for a 



related and quite detailed discussion, see 



35, 



36J). Nevertheless, if in Eq. ([65]) we set C2 = 0, then 



the results coming from both schemes are similar, at least in the sense that, for low curvatures, 
they behave as constant, while for large curvatures the behavior is in both cases proportional to 
R 2 . 

For the novel case of the reconstruction of a YM theory, we followed the same procedure, by 
considering also the same example for both reconstruction schemes. The first one yields the result 



A 1 



1 



In 



hi(l-hi) 
a 2 t 2 



(66) 



while the second scheme gave the following one 



(/H-l)(l-frl) 

3 

1 



5 ln \ — 



(67) 



I In 



MMi) 

cPti 



As before, in the new situation considered in this paper of reconstructing a YM theory, it also 
happens that action (|47j) expresses a more extensive class of theories than the action given by 
(|3ip . and it is again for this reason that more solutions are obtained for the scheme based on an 
auxiliary scalar field. Moreover, with the help of this example we could see explicitly that there is 
in fact a very interesting coincidence between the result obtained in (|66p and one of the results of 
(|67p . This finding here further supports the point of view that the Einsteinian and the Jordanian 
frame descriptions actually lead to two physically different theories, making thus clear the physical 
non-equivalence of the two frames as discussed in the first and in the third references of [29(. In 
view of the strong and still on going discussion about this issue in the specialized literature, this 
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additional piece of evidence is very valuable. Even much more because it comes from a theory that 
it is way closer to physics than the ones considered previously. 

Also important is to remark that, sometimes, it is actually more convenient to use one scheme 
instead of the other, because the final result — after following the way to find the modified or non- 
minimal gravity that reproduces a given desired cosmology — may be definitely easier to obtain and 
to interpret in one of the two schemes. To repeat, although these conclusions may not seem really 
new, since they were already derived in more simplified situations, it is actually comforting (and 
rather non-trivial) to see that they continue to be valid in much more realistic situations, from 
the point of view of physics, as corresponding to the actions here considered, involving Yang- Mills 
fields. 
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